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Abstract. This semi-expository paper surveys results concerning three classes of 
orthogonal polynomials: in one non-hermitian variable, in several isometric non- 
commuting variables, and in several hermitian non-commuting variables. The em- 
phasis is on some dilation theoretic techniques that are also described in some details. 



1. Introduction 

In this semi-expository paper we deal with a few classes of orthogonal polynomials 
associated to polynomial relations on several noncommuting variables. Our initial 
interest in this subject was motivated by the need of more examples related to [THj . 
On the other hand there are transparent connections with interpolation problems in 
several variables as well as with the modeling of various classes of nonstationary systems 
(see for a list of recent references), which guided our choice of topics. Thus we do 
not relate to more traditional studies on orthogonal polynomials of several variables 
associated to a finite reflection group on an Euclidean space or other types of special 
functions of several variables, for which a recent presentation could be found in [21], 
instead we are more focused on results connected with various dilation theoretic aspects 
and Szego kernels. 

Our aim is to give an introduction to this point of view. We begin our presentation 
with a familiar setting for algebras given by polynomial defining relations and then we 
introduce families of orthonormal polynomials associated to some positive functionals 
on these algebras. Section 3 contains a discussion of the first class of orthogonal 
polynomials considered in this paper, namely polynomials in one variable on which 
there is no relation. This choice is motivated mainly by the fact that we have an 
opportunity to introduce some of the basic dilation theoretic techniques that we are 
using. First, we discuss (in a particular case that is sufficient for our goals) the structure 
of positive definite kernels and their triangular factorization. Then these results are 
used to obtain recurrence relations for orthogonal polynomials in one variable with no 
additional relations, as well as asymptotic properties of these polynomials. All of these 
extend well-known results of Szego. We conclude this section with the introduction of 
a Szego type kernel which appears to be relevant to our setting. 

In Section 4 we discuss orthogonal polynomials of several isometric variables. Most of 
the results are just particular cases of the corresponding results discussed in Section 3, 
but there is an interesting new point about the Szego kernel that appears in the proof of 
Theorem 4.1. We also use a certain explicit structure of the Kolmogorov decomposition 
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of a positive definite kernel on the set of non-negative integers in order to produce 
examples of families of operators satisfying Cuntz-Tocplitz and Cuntz relations. 

The final section contains a discussion of orthogonal polynomials of several non- 
commuting hermitian variables. This time, some of the techniques described in Sec- 
tion 3 are not so relevant and instead we obtain three-terms recursions in the traditional 
way, and we introduce families of Jacobi matrices associated to these recursions. Many 
of these results can be proved by adapting the classical proofs from the one scalar 
variable case. However, much of the classical function theory is no longer available 
so we present some proofs illustrating how classical techniques have to be changed or 
replaced. Also some results are not presented in the most general form in the hope 
that the consequent simplifications in notation would make the paper more readable. 



2. Orthogonal polynomials associated to polynomial relations 

In this section we introduce some classes of orthogonal polynomials in several vari- 
ables. We begin with the algebra Vn of polynomials in noncommuting variables 
Xi,...,Xjv with complex coefficients. Let ¥^ be the unital free semigroup on 
generators 1, . . . , A^. The empty word is the identity element and the length of the 
word (7 is denoted by |(7|. It is convenient to use the notation X^r — Xi^ . . . Xi^ for 
cr = ii . . .ifc G F^. Thus, each element P G Vn can be uniquely written in the form 
P = J2ae¥+ '^-crXfj, with Co- 7^ for finitely many cr's. 

We notice that Vn is isomorphic with the tensor algebra over C^. Let (C^)'^'^ denote 
the /c-fold tensor product of with itself. The tensor algebra over is defined by 
the algebraic direct sum 

If {ci, . . . , Bat} is the standard basis of C^, then the set 

{1} U {{Ci, ® . . . ® J 1 < ii, . . . , 4 < AT, /c > 1} 

is a basis of T(C^). For cr = ii . . . we write e„ instead of e^j ® . . . ® e.j^, and the 
mapping X^r cf E F^, extends to an isomorphism from Vn to T(C^), hence 

Vn ^ T(C^) 

It is useful to introduce a natural involution on V2N as follows: 

X^ — XN+k, k — 1, . . . , N, 
X+^Xi_N, l = N+l,...,2N; 

on monomials, 

i^h ■ ■ ■^ik)'^ = ■ ■ ■ 

and finally, ii Q = Z]aeF+^ c^X^, then = SaeF+y CaX^. Thus, V2N is a unital, 
associative, *-algebra over C. 

We say that A C V2N is symmetric if P e ^ implies cP"*" e A for some c e C — {0}. 
Then the quotient of V2N by the two-sided ideal generated by A is an associative 



algebra 71{A). Letting vr = vr^ : V2N 71{A) denote the quotient map then the 
formula 

(2.1) 7r(P)+ = 7r(P+) 

gives a well-defined involution on 71{A). Usually the elements of A are called the 
defining relations of the algebra 71{A). For instance, 7^(0) = 'P2N, 

n{{Xk-x^ \ k = i,...,N}) = VN, 

n{{XkXi - XiXk I A;, / = 1, . . . , 2A^}) ^ ^(C^^), the symmetric algebra over C^^, 

n{{XkXi + XiXk I A;, / = 1, . . . , 2A^}) ^ A(C2^), the exterior algebra over C^^. 
Examples abound in the literature (for instance, see [20] , lU]; EH]). 

There are many well-known difficulties in the study of orthogonal polynomials in 
several variables. The first one concerns the choice of an ordering of In this paper 
we consider only the lexicographic order ^, but due to the canonical grading of 
it is possible to develop a basis free approach to orthogonal polynomials. In the case 
of orthogonal polynomials on several commuting variables this is presented in |21j . 
A second difficulty concerns the choice of the moments. In this paper we adopt the 
following terminology. A linear functional (f) on TZ{A) is called q-positive (q comes from 
quarter) if 0(7r(P)+7r(P)) > for all P EVn- In this case, 0(7r(P)+) = 0(7r(P)) for 
P eVn and 

|0(7r(Pi)+7r(P2))|2 < </)(7r(Pi)+7r(Pi))0(7r(P2)+7r(P2)) 
for Pi, P2 G Vn- Next we introduce 

(2.2) (7r(Pi), 7r(P2))^ = 0(7r(P2)+7r(Pi)), Pi, P2 G Vn- 

By factoring out the subspace A/^ = {vr(P) | P G Vn, (7r(P), 7r(P))0 = 0} and com- 
pleting this quotient with respect to the norm induced by ()2.2j) we obtain a Hilbert 
space H^. 

The index set G{A) C of A, is chosen as follows: if a G G{A), choose the next 
element in G{A) to be the least /? G F^ with the property that the elements 7i{Xa'), 
a' :< a, and 7r(X^) are linearly independent. We will avoid the degenerate situation in 
which 7r(l) = 0; if we do so, then G G{A). Define = TT{Xa) for a G G{A). For 
instance, G(0) = F^, in which case = X^, a G F^. Also, 

G{{XkXi - XiXk I fc, / = 1, . . . , 2N}) = {^1 . . . 2fc G F+ I ^1 < . . . < ^fc, > 0}, 

and 

G{{XuXi + XiXu I A;, / = 1, . . . , 2N}) = {zi . . . G F+ | ^i < . . . < z,, < A; < A^} 

(we use the convention that for A: = 0, ii . . . ik is the empty word). 
We consider the moments of to be the numbers 

(2.3) s^^p = (p{F^Fp) = {Fp,F^)^, a,(3e G{A). 

The kernel of moments is given by K^{a,l3) = Sa,j3, a,(3 E G{A). We notice that 
is a q-positive functional on 1Z{A) if and only if is a positive definite kernel on 
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G{A). However, does not determine uniquely. One typical situation when 
determines is {X^ — \ k = 1, . . . , N} G A; a more general example is provided 
by the Wick polynomials, 

N 

XkX^ — ak,i5k,i — ^ c^^l^X^Xn, k,l — 1,. . . ,N, 

m,n=l 

where a^^i, c^f are complex numbers and 5k,i is the Kroncckcr symbol. 

The moment problem is trivial in this framework since it is obvious that the numbers 
Sa,f3i a,l3 E G{A), are the moments of a q-positive functional on 71{A) if and only if 
the kernel K(a,P) = Sa,p, q;,/5 G G{A), is positive definite. 

We now introduce orthogonal polynomials in TZ{A). Assume that (p is strictly q- 
positive on Ti{A), that is, 0(7r(P)+7r(P)) > for 7r(P) 7^ 0. In this case = {0} and 
'^{'Pn) can be viewed as a subspace of Tiff,. Also, {Fa}aeG{A) is a linearly independent 
family in Ti,^ and the Gram-Schmidt procedure gives a family {<Pa}aeG{A) of elements 
in 71 {Vn) such that 

(2.4) ipa = ^aa,pFf^, aa,a > 0; 

(2.5) {Va,Vi3)ct>^^a,p, a,peG{A). 

The elements (fa, ex & G{A), will be called the orthonormal polynomials associated to 
(f). An explicit formula for the orthonormal polynomials can be obtained in the same 
manner as in the classical, one scalar variable case. Thus, set 

(2.6) Da = det [sa',p']a',p'<a > 0' « ^ ^(-4)' 

and from now on r — 1 denotes the predecessor of r with respect to the lexicographic 
order on F]^, while cr + 1 denotes the successor of a. We have: </?0 = •^00^^ ^iid foi" 
^ a, 

F^ ... Fa \ 

with an obvious interpretation of the determinant. In the following sections we will 
discuss in more details orthonormal polynomials associated to some simple defining 
relations. 

3. No RELATION IN ONE VARIABLE 

This simple case allows us to illustrate some general techniques that can be used in 
the study of orthonormal polynomials. We have A = ^ and N = 1, so TZ{A) = Vi. The 
index set is No, the set of nonnegative integers, and F„ = , n G Nq. The moment 
kernel of a q-positive functional on Vi is K^{n, m) — 0((X{*)+X{"), n,m E No, and we 
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(2.7) 



1 



V Da-lDa 



det 



notice that there is no restriction on other than being positive definite. We now 
discuss some tools that can be used in this situation. 



3.1. Positive definite kernels on Nq. We discuss a certain structure (and parametriza- 
tion) of positive definite kernels on Nq. The nature of this structure is revealed by 
looking at the simplest examples. First, we consider a strictly positive matrix 

1 a 



a 1 



a e 



This matrix gives a new inner product on by the formula 

{x,y)s ^ {Sx,y), a;,^/GM^ 

where (■, ■) denotes the Euclidean inner product on R^. Let {61,62} be the standard 
basis of M^. By renorming with (•, ■)s the angle between ei and 62 was modified to 
the new angle 9 — 0{ei, 62) such that 

(61,62)5 



(3.1) 



cos 61(61, 62) = 



|ei||5||e2||5 

We can visualize the renormahzation process by giving a map T5 : R^ — > with the 
property that {Tsx,Tsy) — {x,y)s for x,y e R^, and it is easily seen that we can 
choose 

1 cos^ 
sine 



Ts 



Co 



Ts 



— >- 






ei 



ei 



Figure 1. Renormahzation in R^ 

We can also notice that Tgei — si and 7562 — f2 — J{cos9)ei, where J {cos 9) is the 
Julia operator, 

cos 9 sin 9 
sin 9 — cos 9 

which is the composition of a reflection about the x-axis followed by the counterclock- 
wise rotation R0 through angle 9. We deduce that 

a — cos9 = (61, /2) = (61, J(cose)6i) = (J(cose)6i, 61). 



J{cos9) 



The discussion extends naturally to the 3x3 case. Thus let 

lab 
S — ale , a,b,c E 
b c 1 



be a strictly positive matrix. A new inner product is induced by on M^, 

{x,y)s = {Sx,y), x,?/eM^ 

and let {ei, 62, 63} be the standard basis of M^. With respect to this new inner product 
the vectors ei, 62, 63 still belong to the unit sphere, but they are no longer orthogonal. 
Thus, 

a = cos 9(61,62) = cos^^i2, 

C = cos 6^(62, 63) = COS 6^23, 

and 

b = cos 9{6i, 63) = cos 6^13. 

This time, the law of cosines in spherical geometry gives a relation between the numbers 
a, b, and c, 

(3.2) b = cos 6*13 = cos 6^12 cos 6^23 + sin 6^12 sin 6*23 cos 9, 

where 9 is the dihedral angle formed by the planes generated by ei, 62 and, respec- 
tively, 62, 63 (see, for instance, [lEl)- Thus, the number b belongs to a disk of center 
cos 6'i2 cos ^23 and radious sin ^12 sin ^23- Once again the renormalization can be visu- 
alized by a map T5 : — > such that {Tsx,Tsy) = {x,y)s- In this case we can 
choose 

1 cos 9i2 cos 9i2 cos 6'23 + sin 9i2 sin 6*23 cos 9 
Ts = sin 6*12 sin 9i2 cos ^23 — cos 6*12 sin 6'23 cos 9 
sin 6'23 sin 9 

and we see that T^ei = ei, Ts62 = f2 = {J{cos9i^2) © l)ei, and 

^563 = /a = iJicos9u) © 1)(1 © J(cose))(J(cos^23) © l)ei. 

In particular, 

(3.3) b = cos 013 = ((J(cos0i,2) © 1)(1 © J {cos 9)) {J {cos 92,3) © l)ei,ei), 
which can be viewed as a dilation formula. 

Now both ()3.2|) and (j3.3|) extend to a strictly q-positive n x n matrix and provide 
a parametrization and therefore a structure for the positive definite kernels on Nq (for 
general results and details see [131, IIHI)- We apply this result to a kernel associated 
to a strictly q-positive functional and obtain that is uniquely determined by 
a family {7fc,j}o<fc<j of complex numbers with the property that |7fcjl < 1 for all 
< k < j. Define dkj = { 

(3.4) 

X (j — A; -|- 1) unitary matrix defined recursively by: Uk^k 



1 — |7fc,jP)^^^- The extension of ()3.3|) mentioned above gives 
Sk,j = s]^ls]'J{Ukj6i, ei), k < j. 



A; + ll 



where Ukj is a {j 
and for k < j, 

(3.5) Uk,j = {J{lk,k+l) © ln~l)(l © J{lk,k+2) © ln-2) . . . (ln-1 © J{-fk,j)){Uk+l,j 

Also J{'ji,m) is the Julia operator associated to 7;,^, that is, 

dl^m 7z,m 
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!)• 



J{li. 



and Im denotes the m x m identity matrix. For instance, we deduce from ()H.4|1 that: 

1/2 1/2 r -I 1 
^01 = So'o s/i [ 1 J 



7oi 






' 1 " 




-7oi . 








so2 = sli^sli^ [ 1 ] 



X 



and 



■S03 



1/2 1/2 
"'OO *33 



X 



X 



7oi 


1 

"01 







"01 


-7oi 
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1 
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do2 
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In particular, we deduce that Sqi = ■Sqq^s]^(^7oi. The next formula is the complex version 
of So2 = -500^322^ (7oi7i2 + c?oi7o2C?i2) • Then, 

503 = Soo^Sg^^ (70l7l2723 + 70lC?127l3C?23 + '^0l702'^12723 
-C?0l7027l27l3C?23 + 'ioiC?02703<^130?23) • 

Explicit formulae of this type can be obtained for each Sk 



allowing to calculate 7^,- from the kernel of moments. 



J' 
see 



as well as inverse algorithms 
[I6J for details. 

A natural combinatorial question would be to calculate the number N{skj) of ad- 
ditive terms in the expression of s^j- We give here some details since the calculation 
of N{sk,j) involves another useful interpretation of formula ()3.4|) . We notice that for 
k>0, ' 

N{sk,k+i) = 1, 
Nisk,k+2) = 2, 
Nisk,k+3) = 5. 

The general formula is given by the following result. 
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N{So2) 
N{S03) 



Theorem 3.1. N{sk,k+i) is given by the Catalan number Ci 



l + l 



21 
I 



Proof. The first step of the proof considers the realization of Skj through a time varying 
transmission hne (or lattice). For illustration we consider the case of sqs in Figure 2. 



A 



B 





Figure 2. Lattice representation for sqs 

Each box in Figure 2 represents the action of a Julia operator, and we see on this 
figure that the number of additive terms in the formula of sqs is given by the number of 
paths from A to B. In it clear that to each path from A to -B in Figure 2 it corresponds 
a Catalan path from C to in Figure 3, that is, a path that never steps below the 
diagonal and goes only to the right or downward. 



€) O O 



o © — e — © 

o o o () 
o o o © 



D 



Figure 3. A Catalan path from C to D 

Thus, each box in Figure 2 corresponds to a point strictly above the diagonal in 
Figure 3. Once this one-to-one correspondence is established, we can use the well- 
known fact that the number of Catalan paths like the one in Figure 3 is given exactly 
by the Catalan numbers. □ 

3.2. Spectral factorization. The classical theory of orthogonal polynomials is intimately 
related to the so-called spectral factorization. Its prototype would be the Fejer-Riesz 
factorization of a positive trigonometric polynomial P in the form P = where 
Q is a polynomial with no zeros in the unit disk. This is generalized by Szego to the 
Szego class of those measures on the unit circle T with logyu' G L^, and very general 
results along this line can be found in |25 and |33j . 

Here we briefly review the spectral factorization of positive definite kernels on the 
set No described in ^3] . For two positive definite kernels Ki and K2 we write Ki < K2 
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if K2 — -ft'i is a positive definite kernel. Consider a family JF = {jF„}„>o of at most 
one-dimensional vector spaces and call lower triangular array a family = {0fcj}fc,j>o 
of complex numbers 6^^- with the following two properties: 0^^- = for k < j and 
each column Cj(6) = [Qk,j]k>o, j > 0, belongs to the Hilbert space ®k>j^k- Denote by 
T-[q{J^) the set of all lower triangular arrays as above. An element of ?ig(jF) is called 
outer if the set {cj(9) | j > A;} is total in ®j>k^j for each > 0. 

It is easily seen that if is an outer triangular array, then the formula 

i^0(fc,j) = cfe(e)*c,(e) 

gives a positive definite kernel on Nq. The following result extends the above mentioned 
Szego factorization and at the same time it contains the Cholesky factorization of 
positive matrices. 

Theorem 3.2. Let K he an positive definite kernel on Nq. Then there exists a family 
^ = {^n}n>o of at most one- dimensional vector spaces and an outer triangular array 
e G nliJ^) such that 

(1) Ke < K. 

(2) For any other family T' = {^n}n>o of at most one- dimensional vector spaces 
and any outer triangular array 0' G Ti-Q^S, T') such that Kq/ < K , we have Kqi < Kq. 

(3) is uniquely determined by (a) and (6) up to a left unitary diagonal factor. 

It follows from (3) above that the spectral factor O can be uniquely determined 
by the condition that Qn,n > for all n > 0. We say that the kernel K belongs to 
the Szego class if inf„>o 0n,n > 0. If {jkj} are the parameters of K introduced in 
Subsection 3.1 then it follows that the kernel K belongs to the Szego class if and only 
if 

(3.6) inf4/,'n^'^.">0- 

n>k 

This implies that J-'n = C for all n > (for details see or jl6j). 

3.3. Recurrence relations. Formula (j2.7|) is not very useful in calculations involving 
the orthogonal polynomials. Instead there are used recurrence formulae. In our case, 
^ = and = 1, we consider the moment kernel K^p of a strictly q-positive functional 
on Vi and also, the parameters {'Jkj} of as in Subsection 3.1. It can be shown that 
the orthonormal polynomials associated to (p obey the following recurrence relations 

(3.7) MXi, I) = ^{X,, I) = s"//', / G No, 
and for n > 1, Z G No, 

(3.8) ifniX,, I) = (Xiy.„_i(Xi, / + 1) - ^i,n+wl_,{Xi, I)) , 

(3.9) ^liX.J) = (-7,,„+,Xiy,„_i(Xi,/ + 1) + ^t.iX^J)) , 
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where ^n{Xi) = </^„(Xi,0) and y^KXi) = ipi{Xi,0). 

Somewhat similar polynomials are considered in , but the form of the recurrence 
relations as above is noticed in [T7j. It should be mentioned that {v5n(-^i, 0}n>o is 
the family of orthonormal polynomials associated to a q-positive functional on Vi with 
moment kernel K\a,P) = Sa+i,f3+i, G Nq. Also, the above recurrence relations 
provide us with a tool to recover the parameters {'Jkj} from the orthonormal polyno- 
mials. 

Theorem 3.3. Let k}^ he the leading coefficient of (pn{Xi,l) . For I G No and n>l, 

1,1+1 ,1+1 
7;,n+; = -<^n(0,0— w T^- 

K,Q . . . fi,^ 

Proof. We reproduce here the proof from [8] in order to illustrate these concepts and 
to introduce one more property of the parameters {'Jkj}- First, we deduce from (|3.8p 
that 

¥^n(0,0 = -y^¥'i-l(0,0, 

while formula ()3.9p gives 

'pUo, I) = ;^vLM /) = ... = .-//^ n 

(J'l,n+l " Cll,p+l 

hence 

n ^ 

<^n(0, /) = -S^l^'^-fi^n+l n 1 • 

p=l 

Now we can use another useful feature of the parameters {'Jkj}, namely the fact that 
they give simple formulae for determinants. Let Dm,i denote the determinant of the 
matrix [sk,j]i<k,j<m- Proposition 1.7 in jlSI, 

m 

(3.10) Di,m = Ylsk,kX n 4p- 

k=l '<i<P<'T* 

One simple application of this formula is that it reveals the equality behind Fisher- 
Hadamard inequality. Thus, for I < n < n' < m, we have 



Dl,n'Dn,m TT ,2 
Di^rn-—^ 11 '^fc.J' 



(fcj)gA 



where A = {{k,j) \ l<k<n<n'<j< m}. Some other applications of ()3.10|) can 
be found in [121, Chapter 8. Returning to our proof we deduce from ()3.10|) that 



,2 _ ^-1 Dl,l+n 
Dl+l,l+n 



p=i 
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(3.11) 7l,n+l = -Vn{0,l). 



Di 



l+n 



l+n 



We can now relate this formula to the leading coefficients k}^. From ()3.8j) we deduce 
that 

n— 1 ^ 

kn = Si^^ i^^\\ - , n > 1, 
and using once again (j3.1(J|) . we deduce 



which concludes the proof. □ 



3.4- Some examples. We consider some examples, especially in order to clarify the con- 
nection with classical orthogonal polynomials. Thus, consider ffist ^ = {1 — X^Xi}. 
In this case the index set is still Nq and if is a linear functional on TZ{A), then the 
kernel of moments is Toeplitz, K^{n + k,m + k) = K^{n, m), m,n,k & Nq. Let be a 
strictly q-positive functional on 71{A) and let {'Jkj} be the parameters associated to K^. 
We deduce that these parameters also satisfy the Toeplitz condition, 'jn+k,m+k = Jn,m, 
n < m, k > 1. Setting 7„ = 'jk,n+k, n > 1, k > 0, and dn = (1 — |7„p)^/^, the recurrence 
relations ()3.8|) . ()3.9p collaps to the classical Szego recursions obeyed by the orthogonal 
polynomials on the unit circle, 

"n+1 

and ^ 

"n+1 

Therefore 7„, n>l are the usual Szego coefficients, 

Another example is given by ^ = {Xi — X^}. In this case the index set is still 
No and the moment kernel of a strictly q-positive functional on TZ{A) will have the 
Hankel property, K^{n., m + k) = K^{n + /c, m), m, n, G Nq. Orthogonal polynomials 
associated to functionals on TZ{A) correspond to orthogonal polynomials on the real 
line. This time, the parameters {7fc,j} associated to moment kernels have no classical 
analogue. Instead there are so-called canonical moments which are used as a counter- 
part of the Szego coefficients (see |2ZI). Also, recurrence relations of type ()3.8p . ()3.9|) 
are replaced by a three term recurrence equation, 

(3.12) Xipn{x) = bn^n+l{x) + a^ipnix) + b„,_iipn-l{x) , 

with initial conditions (p-i = 0, (po = 1 (|32|). Definitely, these objects are more useful 
(for instance, it appears that no simple characterization of those {'Jkj} corresponding to 
Hankel kernels is known). Still, computations involving the parameters {'Jkj} niight be 
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of interest. For instance, we show here how to calculate the parameters for Gegenbauer 
polynomials. For a number A > — |, these are orthogonal polynomials associated 

to the weight function w{x) = B{^,X + — x'^)^~^ on (—1,1) {B denotes the 

beta function). We use the normalization constants from j^, thus the Gegenbauer 
polynomials are 

where (x)„ is the Pochhammer symbol, {x)q = 1 and (x)„ = YYk=ii^ ~^ ^ ~ ^) ^"^^ n>l. 
We have: 

- ^ t (P'(x)? (1 - x'V--2dx - + 

" X + i) ^^^^ ^ ^ - 2(2A + lUn + A) 

and the three term recurrence is: 

(see Ch. 1). We now let {p^{x,0) denote the orthonormal polynomials associated 
to the weight function w, hence ip^(x,0) = -^=P^(x). From the three term relation 

we deduce 



n\{n + 2X) ^J-k-l + 2X 



and also, the leading coefficient of V5„(a;, 0) is 



,A,o_ in + 2X)n 2{2X + lUn + X) 
" 2^{X + l)J n\{n + 2X) ' 

In order to compute the parameters {7^ j} of the weight function w we use Theo- 
rem ESI Therefore we need to calculate the values (/?^(0, 1) and k^, n > 1, I > where 
ky denotes the leading coefficient of (p^{0,l). The main point for these calculations 
is to notice that {ip^{x,l)}n>o is the family of orthonormal polynomials associated to 
the weight function x'^''w{x). These polynomials are also classical objects and they 
can be found for instance in [21] under the name of modified classical polynomials. A 
calculation of the modified Gegenbauer polynomials can be obtained in terms of Jacobi 
polynomials. These are orthogonal polynomials associated to parameters a, /3 > 1 and 
weight function 

2-"-^-^5(a + 1, /5 + 1)^^(1 - a;)"(l + xf 
on (—1, 1) by the formula 

p!C'^\x) = ^(1 - x)-°(l + x)-^ — {l - x)"+"(l + xr^\ 
According to P^, Sect. 1.5.2, we have 

^l{x,l)=C,nPn~'"'~H'^x'-^) 
12 



and 



where c„ is a constant that remains to be determined. But first we can aheady notice 
that the above formulae give f2n+ii^y = O5 so that ■y^2n+i+i — 0- 



Theorem 3.4. Forn,l > 1, 



n 



X+l+k-l 




and 



,X,l_ (A + /)2n /(A + 1) 



.A,I _ (A + /)2„+l / (A+1), 



where 



and 



(A + 1)JA + /UA + /) 



,A. (A + |)JA + /Ui(A + /) 



^! (^ + i)n+i (A + / + 2n + l)- 
Proof. It is more convenient to introduce the polynomials 

and again by classical results that can be found in ||21j, we deduce 

1 

1 



.2 f i^Hu y gg+jA+i) A,; 

V(A+Ouy B(|,A+|) '*2n> 



where 



A._ (A+|)n(A + /)„(A + /) 

!(Z + |)„(A + Z + 2n)- 
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Using that ^^X'^'^X^? = ttttv? we deduce 



V (2); ^2; 

The calculation of ^P2ni^y^) reduces to the calculation of 6*2^^(0) which can be easily 
done due to the three term relation 



Thus we deduce 



^^2^^2(0) = -^^<(0), 



and by iterating this relation and using that Cq' (0) = 1, we get 

- 1 



¥^2^(0,/) 



\n+l 




The leading coefficient of (/52n(^; can be obtained from the corresponding formula in 
1. Thus, 



and 



'^2n+l 




where 



"'2n+l 



(A + l)JA + 0.+i(A + /) 

i} + \)n+i (A + ^ + 2n + l)' 

□ 

Now the parameters {7^^} can be easily calculated by using Theorem 13 .41 Of course, 
the explicit formulae look too complicated to be recorded here. 

3.5. Asymptotic properties. In the classical setting of orthogonal polynomials on the 
unit circle there are several remarkable asymptotic results given by Szego. Let be a 
measure in the Szego class, and let {v?n}n>o be the family of orthonormal polynomials 
associated to /z. Then, the orthonormal polynomials have the following asymptotic 
properties: 

(3.13) y^n ^ 
and 

(3.14) 4 ^ ^z^' 
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where 0^ is the spectral factor of /i and the convergence is uniform on the compact 
subsets of the unit disk D. The second limit (I3.14j) is related to the so-called Szego 
limit theorems concerning the asymptotic behaviour of Toeplitz determinants. Thus, 

det Tn 1 

detT„_i ~ |(^|(0)|2' 

where T„ = [sj_j]"j^Q and {sk}kez is the set of the Fourier coefficients of /i. As a 
consequence of the previous relation and (j3.14|) we deduce Szego's first limit theorem, 

detT 1 f^'^ 

(3.15) hm = |e,(0)p = exp(- / log f,'{t)dt). 

n^oo det ln-1 277 Jq 

The second (strong) Szego limit theorem improves (j3.15j) by showing that 



where g{fi) is the limit in (j3.15|) and a is the planar Lebesgue measure. These two 
limits (j3.15j) and (j3.16|) have an useful interpretation in terms of asymptotics of angles 
in the geometry of a stochastic process associated to fi (see (23]) and many important 
applications. We show how these results can be extended to orthogonal polynomials 
on Vi. The formulae ()3.8p and ()3.9|) suggest that it is more convenient to work in 
a larger algebra. This is related to the so-called Toeplitz embedding, see jTH], [231 • 
Thus, we consider the set £ of lower triangular arrays a = [ak,j]k,j>o with complex 
entries. No boundedness assumption is made on these arrays. The addition in C is 
defined by entry-wise addition and the multiplication is the matrix multiplication: for 
a = [ak,j]k>j, b = [bk,j]k,j>o two elements of C, 

{oh)k,j = y^^akAj^ 

which is well-defined since the sum is finite. Thus, C becomes an associative, unital 
algebra. 

Next we associate the element of C to the polynomials (y9„(Xi, /) = "^2=0 '^n fc^i ' 
77,, Z > 0, by the formula 

k<j- 

similarly, the element of C is associated to the family of polynomials (pl^{Xi,l) = 
Ym:=o bn^k^i, n, / > 0, by the formula 

(3.18) ($«K. = {|;-^- l^i 

We notice that the spectral factor 0^ of is an element of C and we assume that 
0,^ belongs to the Szego class. This implies that is invertible in C for all n > 0. 
Finally, we say that a sequence {a„} C C converges to a e £ if {{an)k,j} converges to 
Qkj for all k,j > (and we write a„ — > a). 
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(3.17) ($„) 



Theorem 3.5. Let belong to the Szego class. Then 



(3.19) $n ^ 

and 

(3.20) ($1)-^ ^ e^. 

We now briefly discuss the geometric setting for the kernel K^. By a classical result 
of Kolmogorov (see [SU]), is the covariance kernel of a stochastic process {fn}n>o C 
L^(/i) for some probability space {X,A4,fi). That is, 

K^{m,n) = / fnlmdfJ'- 
Jx 

We can suppose, without loss of generality, that {/n}n>o is total in L^(/x) and for p < g 
we introduce the subspaces £p^q given by the closure in L'^i^jj) of the linear span of 
{fkYk=p- The operator angle between two spaces £i and £2 of L'^{lj) is defined by 

B{£r,£2)=Ps,Pe,Ps,, 

where Ps^ is the orthogonal projection of L'^{n) onto £1. Also define 

A(^i,^2) = / -5(^1,^2). 

We associate to the process {/n}„>o a family of subspaces Tir^q of -L^(/i) such that 
Tir^q is the closure of the linear space generated hj fk, r < k < q and we consider a 
scale of limits: 

(3.21) S - lim A(7Yo,n,^n+l,r) = A(7^o,n,^n+l,oo) 

r— »oo 

for n > 0, and then we let n ^ 00 and deduce 

(3.22) S— lim A(7io,n, "^^n+l.oo) = A(7io,oo, n„>o7in,oo), 

n— >oo 

where s — lim denotes the strong operatorial limit. 

We then deduce analogues of the Szego limit theorems ()3.15|) and (j3.16|) by expressing 
the above limits of angles in terms of determinants. This is possible due to ()3.10|) . 

Theorem 3.6. Let (p belong to the Szego class. Then 

(3.23) '^''^ = Sr,r det A(?ir,r, 'Hr+l 



P>r+l,q kg-r(0,r)|2 

and 

(3.24) lim = Sr,r det A{Hr,r,T-tr+l,oo) = |60(r,r)|^ = Sr.r J]^ C?r,r+i- 

J>1 

// we denote the above limit by gr and 



L = lim TT dl > 0, 



n— >oo 

0<fc<n<j 
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then 



(3.25) 



lim 



Do,n _ 1 _ 1 



■oo 



nr=o 9i det A (7^0,00, n„>oH n,oo / 



Details of the proofs can be found in |HI. 

3.6. Szego kernels. The classical theory of orthogonal polynomials is intimately related 
to some classes of analytic functions. Much of this interplay is realized by Szego kernels. 
Here we expand this idea by providing a Szego type kernel for (a slight modification 
of) the space TiKJ-") which is viewed as an analogue of the Hardy class on the unit 
disk. We mention that another version of this idea was developed in j2] (and recently 
applied to a setting of stochastic processes indexed by vertices of homogeneous trees 
in [S]); see also 0, 0- The difference is that TCq{J-') is larger than the space U2 of [2] 
and also, the Szego kernel that we consider is positive definite. 

We return now to the space Hl{J^) introduced in Subsection 3.2. Its definition 
involves a family of Hilbertian conditions therefore its natural structure should be that 
of a Hilbert module (we use the terminology of |2H])- Assume J-'„ = C for all n > and 
consider the C*-algebra V of bounded diagonal operators on ©„.>o-^n = ^^(Nq). For a 
sequence {dn}n>o we use the notation 



so that diag {{dn}n>o) belongs to V if and only if sup„>Q \dn\ < 00. We now define the 
vector space 



turning into a Hilbert P-module. As a Banach space with norm ||0|| = 

11(6, 6)11^/^ the space H^{T) coincides with /°°(No, /^(Nq)), the Banach space of bounded 
sequences of elements in /^(Nq). A similar construction is used in |TT] for a setting of 
orthogonalization with invertible squares. 

Next we introduce a Szego kernel for 7i^(jF). Consider the set 



do 
di 



1 







diag ({c/„}„>o) 








Bl = {{Zn}n>0 C C I sup \Zn\ < 1} 



n>0 
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and for z — {2;n}n>o G -Bi, © G ^{^{T), notice that 



Q{z) = diag I {Qn,n + ^ Qk,nZk-l ■ ■ ■ Zn}n>0 1 
\ k>n / 



is a well-defined element of V. Also, for z e Si, we define 

1 
^0 1 

ZqZi Zi 1 
ZqZiZ2 Z1Z2 Z2 



which is an element of and the Szego kernel is defined on Bi by the formula: 

(3.26) S{z,w) = {Syj,Sz), z,w e Bi. 

Theorem 3.7. S is a positive definite kernel on Bi with the properties: 

(1) e{z)^{0,s,), een\T),zeB,. 

(2) The set {S,D \ zeBi,D eV} is total in rL^{T). 

Proof. Take zi, . . .Zm £ Bi and after reshuffling the matrix [S{zj, zi)]^^^^ can be written 
in the form 

Since each matrix [c„(S'^^)*Cn(S'^;J]'"^_^ 



Cn ('S'zi) 


* 








. Cn{Sz^) _ 



is positive, we 



conclude that 5" is a positive definite kernel on Bi. Property (1) oi S follows directly 
from definitions. For (2) we use approximation theory in L°° spaces in order to reduce 
the proof to the following statement: for every element © e T-C^{J-') for which there 
exists A C N such that c„(0) = h e /^(No) for n G A and c„(0) = for n ^ A, 
and for every e > 0, there exists a linear combination L of elements SzD. z G Bi, 
D E V, such that sup„>o l|cn(0 ^ -^)|| < This can be achieved as follows. Since 
the set {(j)w{z) = I ^ -^i} is total in the Hardy space H"^ on the unit disk, we 
deduce that there exist complex numbers Ci, . . ., and Wi, . . ., w^, \wk\ < 1 for all 
k — 1, . . .m, such that 

1 



\\h 



k=l 



Wk 

wi 



< e. 



Then define Zk — {■u;jk,„}„>o ior k — 1, . . . ,m, where Wk,n — Wk for all n > 0. So 
Zk E Bi. Also define — 1 ior n E A and d^ — ior n ^ A, and consider 

m 

L = '^CkSz^di&g{{d^}) . 



k=l 
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We deduce that L G H'^{J^), Cn{L) = for n ^ A, and Cn{L) = Yl'k=i^k 
n G A, so that 

SUP„>0 ||Cn(0 - ^)} = SUp„>o ||c„(e) - C„(L)} 

= max{sup„g^ \\h - c„(L)||, sup„^^ l|c„(^)||} 
1 " 

□ 

4. Several isometric variables 

In this section we discuss orthogonal polynomials in several variables satisfying the 
isometric relations X+X^ = 1, = 1, . . . , X. We set ^ = {1 - X+X^ | A; = 1, . . . , X} 
and notice that the index set of A is F^. Also if is a linear functional on TZ{A) then 
its kernel of moments is invariant under the action of on itself by juxtaposition, 
that is, 

(4.1) K^{Ta, ra') = K^a, a'), r, a, a' G F+ . 

In fact, a kernel K satisfies (14. Ij) if and only if i^' = for some linear functional 
on TZ{A). Positive definite kernels satisfying (j4.H) have been already studied, see for 
instance jH] and references therein. In particular, the class of isotropic processes on 
homogeneous trees give rise to positive definite kernels for which a theory of orthogonal 
polynomials (Levinson recursions) was developed in jlOj. Here we discuss in more 
details another class of kernels satisfying ()4.1|) which was considered, for instance, in 

m- 

4.1. Cuntz-Toeplitz relations. Consider the class of positive definite kernels with prop- 
erty ()4.H) and such that 

(4.2) K{cr, r) = if there is no a G F^ such that cr = ra or r = aa. 

We showed in that K has properties (j4.1|) and (j4.2j) if and only if = for 
some q-positive functional on 71{Act), where Act = {1 — X^X^ \ k = 1, . . . , N} U 
{X^X;, k,l = 1, . . . ,N,k 7^ I}. The relations in Act are defining the Cuntz-Toeplitz 
algebra (see [221 details). The property ()4.2|) shows that K is quite sparse, therefore 
it is expected to be easy to analyse such a kernel. Still, there are some interesting 
aspects related to this class of kernels, some of which we discuss here. 
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Wk 



for 



Let be a strictly q-positive kernel on 71{Act) and let be the associated kernel 
of moments. Since the index set of Act is still F^, and this is totally ordered by the 
lexicographic order, we can use the results described in Subsection 3.1 and associate to 
a family {■y^.^}^^^ of complex numbers with |7o-,r| < 1, uniquely determining by 
relations of type (j3.4p . It was noticed in [T7] that has properties (|4.1|) and (|4.2p if 
and only if jTa,Ta' = 1a,a' and 7o-,r = if there is no a G such that a = ra or r = aa. 
The main consequence of these relations is that K^f, is uniquely determined by ■ja- = 70,0-, 
a e F]^ — {0}. We define = {1 — 170-^)^''^. The orthogonal polynomials associated 
to (p satisfy the following recurrence relations (see [T^ for details): = = s^^^'^ 
and for /c G {1, . . . , A^}, a G F^, 

1 u 

(4.3) ^ka = —{XkV^ - 7fca<^L-l)' 

"fccr 

(4.4) i = ;^^(-7fc.^fc^. + V^L^i)- 

The results corresponding to Theorem 13.51 and Theorem 13.61 can be easily obtained 
(see jH]), but the constructions around the Szego kernel are more interesting in this 
situation. Thus, there is only one Hilbertian condition involved in the definition of 
1-Lq{T) in this case. In fact, it is easy to see that Ho{J^) can be identified with the 
full Fock space /^(F^), the space over F^. Now, concerning evaluation of elements 
of 7^0 (-^), if we are going to be consistent with the point of view that the "points 
for evaluation" come from the unital homomorphisms of the polynomial algebra inside 
TioiJ^), then we have to consider an infinite dimensional Hilbert space S and the set 

N 

B,{£) = {Z= (Zi, ...,Z!,)e C{£f I J2 ZkZl < /}. 

k=l 

For a = ii . . Ak G F^ we write Z(j instead of Zi^ . . . Zi^. Then we define for G 
/2(F+) ®^ and Z G Bi{£), 

which is an element of the set C{S) of bounded linear operators on the Hilbert space 
S. Next, for Z G Bi{£) we define Sz : S ^ /2(F+) ® f by the formula: 

Szf = 5^ ® (^.)7, / e 

Then Sz G C{S,l'^(¥j^) (g) S) and we can finally introduce the Szego kernel on Bi{£) 
by the formula: 

S{Z, W) = S*zSw, Z,W e Bi{£). 
Theorem 4.1. S is a positive definite kernel on Bi{S) with the properties: 

(1) eiz) = S*ze, Q ef{¥+)®£,Z eBi{£). 

(2) The set {Szf \ Z e Bi{8)J & 8} is total m /2(F+) ® £. 

20 



Proof. The fact that 5* is positive definite and (1) are immediate. More interesting is 
(2) and we reproduce here the proof given in [T7j. Let / = {fa}o-£¥+ element 
of (g) S orthogonal to the linear space generated by {Szf \ Z G Bi{£), f e S}. 

Taking Z = 0, we deduce that = 0. Next, we claim that for each a E — {0} 
there exist 

Zi = iZ[,...,Z'^)eB,i£), / = !,.. .,2|a|, 



such that 
and 



range 



7*1 r7*2\u\ 

'„ . . . Z/o- 



Z| = for all Ty^a, |r| > / = 1, . . . , 2|a|. 

Once this claim is proved, a simple inductive argument gives / = 0. In order to 
prove the claim we need the following construction. Let {^^jYlj^x be the matrix units 
of the algebra M„ of n x n matrices. Each e^- is an n x n matrix consisting of 1 in 
the {i,j)th entry and zeros elsewhere. For a Hilbert space £i we define E^j = e^- le^ 

and we notice that E^jEJ^i = SjkE^^ and E*^ = E^j. Let a = ii . . .ik so that S = f:®^''^' 
for some Hilbert space Si (here we use in an essential way the assumption that S is of 
infinite dimension). Also, for s = 1, . . . , A^, we define = {/ G {1, . . . , A;} | ik+i-i = s} 
and 

1 



^7 — ^ Er+!)-l,r+pJ S — 1, . . . , iV, p — 1 



reJs 

We can show that for each p G {1, . . . , |(t|} 



(4-5) Z*f = -^Ef;[^, 



2k 

(4.6) ZP = for r^a, |r| > |(t|. 

We deduce 

EN yp y*p _ 1\^N p2|o-| Tp2\a\ 

s=l 2 2^s=l l^r&Js -^r+p,r+p-l-^r+p-l,r+p 

2 Z^s=l Z^r£Js ^r+p,r+p 

2 Z—ir=l r+p,r+p ^ ! 

hence Z^ G Bi{£) for each p = 1, . . . , |(t|. For each word t = ji . . . jfc ^ ^ti " {0} 
deduce by induction that 



(4.7) Z*" ...Z*" = E'I'I 



^ / J r+p— l,r+p+A:-l' 



where A, = n^J( J,,_, - p) C {1, . . . , A^} and Jj,_^ - p = {/ - p | / G J.,_ J. 

We show that = {1} and A,- = for t ^ a. Let g G A^. Therefore, for any 
p G {0, . . . , — 1} we must have q + p E Jj^-p o^' ik+i-q-p = jk-p- For p = — 1 we 
deduce ji = i2-q and since 2— g > 1, it follows that g < 1. Also g > 1, therefore the only 
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element that can be in A,- is g = 1, in which case we must have r = a. Since / G Ji^j^^_i 
for each Z = 1, . . . , /c — 1, hence A„ = {1} and A,- = for r 7^ a. Formula (j4.7|) implies 
()4.5|) . In a similar manner we can construct a family Z^, p = |cr| + 1,...,2|(t|, such 
that 

1 



zy = 



and 



Thus, for s = 1 N, we define 



+p—k,r+p—k—ly 



2k 

a, 

{I G {1, 
s = 1, 



2\a\ 
p+k,pi 



\t\ > 










P 



N, p= |ct| + l,...,2|a|. 



-eKs 



Now, 



7*1 



Z*2|ct 
o- 



E- 



2|ct| 
l,fc+l 



whose range is £. This concludes the proof. 



rp2\a\ rp2\a\ 
-^k,2k -^k+1,1 



E. 



2\a\ 
2k,k 



□ 



It is worth noticing that property (2) of 5* is no longer true if S is finite dimensional. 
In fact, for £ of dimension one the set {Szf \ Z G Bi{£),f G £} is total in the 
symmetric Fock space of (see 0). 



4-2. Kolmogorov decompositions and Cuntz relations. This is a short detour from 
orthogonal poljTiomials, in order to show a construction of bounded operators satisfying 
the Cuntz- Toeplitz and Cuntz relations, based on parameters {7cr}a-gF+-{0} associated 
to a positive definite kernel with properties (j4.H) and (j4.2|) . 

First we deal with the Kolmogorov decomposition of a positive definite kernel. This is 
a more abstract version of the result of Kolmogorov already alluded to in Subsection 3.5. 
For a presentation of the general result and some applications, see |221, ISO]- Here we 
consider : Nq x Nq — C a positive definite kernel and let {7fc,j} be the family of 
parameters associated to K as in Subsection 3.1. In addition we assume K{j,j) = 1 
for j > 0. This is not a real loss of generality and it simplifies some calculations. We 
also assume I'jkjl < 1 for all k < j. Then we introduce for < /c < j the operator Vkj 
on /^(No) defined by the formula 

Vk, = {J{lk © l„_i)(l © J{lk,k+2) © ln-2) • • • (ln-1 © Jilkj)) © 

and we notice that 

(4.8) Wk = s- lim Vk,j 

is a well-defined isometric operator on Z^(No) for every A; > 0. If we define ^(0) = J/C 
and V{k) = WqWi . . . Wk-i/C for k > 1, then we obtain the following result from pTS] . 
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Theorem 4.2. The map V" : No — /^(Nq) is the Kolmogorov decomposition of the 
kernel K j in the sense that 

(1) K{3,l) = {V{l),V{3)), J,/eNo. 

(2) The set {V{k) | G No} is total m /^(No). 

It is worth noticing that we can write exphcitly the matrix of Wk'- 

7fc,fc+l dk^k+llk,k+2 dk,k+ldk,k+2'yk,k+3 

dk,k+i -7fc,fe+i7fc,fc+2 -lk,k+idk,k+2'lk,k+3 ■ ■ ■ 

dk,k+2 -lk,k+2lk,k+3 

■ dk^k+3 

■ 



We also see that Theorem 3.7 and Theorem 4.1 produce various kinds of Kolmogorov 
decompositions for the corresponding Szego kernels. Based on a remark in jT^j, we 
use Theorem 14. 21 in order to obtain some large families of bounded operators satisfying 
Cuntz-Toeplitz and Cuntz relations. Thus, we begin with a positive definite kernel K 
with properties (j4.H) and (|4.2p . For simplicity we also assume iir(0, 0) = 1 and let 
{7o-}o-gf+-{0} tie the family of corresponding parameters. In order to be in tune with 
the setting of this paper, we assume |7o-| < 1 for all cr e — {0}. Motivated by the 
construction in Theorem 14.21 we denote by CaiWo), o G F^i^ — {0}, the columns of the 
operator Wq. Thus, 



ci(l^o 



7i 
d\ 




and Co- (Wo) 



dx - . . c?o-i7fT 
-7i(i2 . . . d„^xla 



da 





for 1 ^ cr. 



We now define the isometry U{k) on /^(F]l^) by the formula: 

U{k) = [Ckr{W^)],^^+, k = l,...,N. 

Theorem 4.3. (a) The family {Ui, . . . ,Uj\f} satifies the Cuntz-Toeplitz relations: U^Ui = 
6k,il, k,l = l,...,N. 

(6) The family {Ui, . . . ,Un} satifies the Cuntz relations UlUi = 5k,il , k,l = 1, . . . , N 
and J2k=i ^kUk — ^ ' ^/ ^^^^2/ ^/ 



(4.9) 



n da=o. 



aGF+-{0} 
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Proof, (a) follows from the fact that Wq is an isometry. In order to prove (b) we need 
a characteriztion of those Wq which are unitary. Using Proposition 1.4.5 in we 
deduce that Wq is unitary if and only if (j4.9j) holds. □ 

It is worth mentioning that if we define V^(0) = //C and V{cr) = U{a)/C, a G 
~ {^}^ where U{a) = U{ii) . . . U{ik) provided that a = ii . . .ik, then V is the 
Kolmogorov decomposition of the kernel K. A second remark here is that the condition 
()4.9p is exactly the opposite of the condition for K being in the Szego class. Indeed, 
it is easily seen that for a positive definite kernel with properties M.lj) and ()4.2|) . the 
condition ()3.6p is equivalent with no-eF+-{0} '^o- > 0. 

5. Several hermitian variables 

The theory corresponding to this case should be an analogue of the theory of or- 
thogonal polynomials on the real line. We set A = {Yk — | /c = 1, . . . , A^} and 
A! = AU {YkYi - YiYk \ k,l = 1, . . . , N}, and notice that 7^(^) = Vn- Also, 7^(^') 
is isomorphic to the symmetric algebra over C^. Orthogonal polynomials associated 
to 71{A'), that is, orthogonal polynomials in several commuting variables were studied 
intensively in recent years, see j2Tj. In this section we analyse the noncommutative 
case. The presentation follows 

Let be a strictly q-positive functional on T]\i{A2) and assume for some simplicity 
that is unital, 0(1) = 1. The index set of A is and let {v5o-}o-gF+ be the orthonormal 
polynomials associated to 0. We notice that for any P,Q E Vn, 

{XkP,Q)^ = 0(Q+XfcP) 

(5.1) = 0(g+x+p) 

which implies that the kernel of moments satisfies the relation Saa,T = ■5o-,/(a)r for 
a, cr, r G ¥j^, where I denotes the involution on given by /(ii . . . ifc) = ik ■ ■ - ii- This 
can be viewed as a Hankel type condition, and we already noticed that even in the 
one dimensional case the parameters {jk.j} of the kernel of moments of a Hankel type 
are more difficult to be used. Therefore, we try to deduce three-terms relations for the 
orthonormal polynomials. A matrix-vector notation already used in the commutative 
case, turns out to be quite useful. Thus, for > 0, we define P„ = [v5o-]|o-|=n ? > 0, 
and P_i = 0. 

Theorem 5.1. There exist matrices An,k (^nd Bn,k such that 

(5.2) XkPn = Pn+lBn,k + PnAn,k + Pn-lP^i,^, k = 1, . . . , N, U > 0. 

Each matrix A„ ^ is a selfadjoint A^" x A^" matrix, while each P„ ^ is an n""^^ x A^" 
matrix such that 

Bn = [ -B„,l . . . Pn,Ar ] 
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is an upper triangular invertible matrix for every n > 0. For n = —1, B^i^k = 0, 
k = 1, . . . , N. The fact that is upper triangular comes from the order that we use 
on F^. The invertibility of 5 is a consequence of the fact that is strictly q-positive 
and appears to be a basic translation of this information. It turns out that there are 
no other restrictions on the matrices ^n,fe, as shown by the following Favard type 
result. 

Theorem 5.2. Let (fa- = '^T<a^(^,T-^r, cr e F^, be elements in Vn such that = 1 
and a^j^cF > 0. Assume that there exists a family {An^k, Bn,k | 't- > 0, k — 1, . . . , N}, of 
matrices such that A* ^ = A^^k cind Bn — [ -B„^i . . . Bn^^ ] is an upper triangular 
invertible matrix for every n>0. Also assume that 
(5.3) 

where [fa]\fj\=-i ~ '^''^^ -B-i,fc = for k = 1,...,N. Then there exists a unique 
strictly positive functional (p on 71{A) such that {^cr}aG¥+ family of orthonormal 

polynomials associated to (p. 

There is a family of Jacobi matrices associated to the three-term relation in the 
following way. For P G 71{A){= Vn), define 

(5.4) ^4>iP)^'^ = P'P'r- 

Since the kernel of moments has the Hankel type structure mentioned above, it follows 
that each ^</,(-P) is a symmetric operator on the Hilbert space Hcf, with dense domain V, 
the linear space generated by the polynomials (pa, E F^. Moreover, for P,Q E Vn, 

*^(pg) = *^(p)*^((5), 

and ^^{P)!) C T>, hence is an unbounded representation of Vn (the GNS repre- 
sentation associated to 0). Also, (p^P) = (\l/(/,(P)l, !)</, for P G Vn- Wc distinguish 
the operators -^k = ^0(>fe), k = 1,...,N, since *^(Ea6F+ c^^^) = EaeiF+ c<7^</,,a, 
where we use the notation ^^^o- = • • • for cr = ii . . . ifc- Let {ei, . . . , Ctv} be the 
standard basis of and define the unitary operator W from P(¥~j^) onto Ti.^ such 
that W{ea-) = (fa, a G F^. We see that W~^T> is the linear space Vq generated by Co-, 
(7 G F^, so that we can define 



Jk^W-''^^,kW, A; = 1,...,7V, 
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on Do- Each Jk is a symmetric operator on "Dq and by Theorem 15. the matrix of (the 
closure of) Jk with respect to the orthonormal basis {eo-}^gp+ is 



Bo^k Aik Blf^ 



-Bi,fe ^2,fc 



We call ( Ji, . . . , Jjv) a Jacobi A^-family on Vq. It is somewhat unexpected that the usual 
conditions on An^k and i?„ ^ insure a joint model of a Jacobi family in the following 
sense. 

Theorem 5.3. Let (Ji, . . . , Jn) a Jacobi N -family on T>q such that A* ^ = An^k one? 
Bn = [ Bn,i ■ ■ ■ Bn,N ] is au uppcT triangular invertible matrix for every n > 0. 
Then there exists a unique strictly q-positive functional on Vn with associated or- 
thonormal polynomials {<^a}a^f+ such that the map 



Proof. First the Favard type Theorem 15 . 21 gives a unique strictly q-positive functional (p 
on Vn such that its orthonormal polynomials satisfy the three-term relation associated 
to the given Jacobi family, and then the GNS construction will produce the required 



One possible application of these families of Jacobi matrices involves some classes of 
random walks on Fj^. Figure 4 illustrates an example for N = 2 and more details are 
planned to be presented in |7|. 



W{e^) =ip,, ae F+, 
extends to a unitary operator from t^{¥^) onto Ti^ and 

Jk = W~^^^^kW, k = l,... 



N. 



W and '^4i^k^ as explained above. 



□ 




Figure 4. Random walks associated to a Jacobi family, N = 2 



We conclude our discussion of orthogonal polynomials on hermitian variables by 
introducing a Szego kernel that should be related to orthogonal polynomials on Vn- 
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Thus, we consider the Siegel upper half-space of a Hilbert space 8 by 

H+{8) = {{Wi ...Wn)^ C{8)'' I W^Wl + . . . + Wn-iW*^.i < ^{Wn - W*^)]. 

We can estabhsh a connection between Bi{8) and H^{8) similar to the well-known 
connection between the unit disk and the upper half plane of the complex plane. Thus, 
we define the Cayley transform by the formula 

C{Z) = ((/ + Z^yZ,, ...,(/ + Z^)-^Zjv_i, 2(/ + Z^)-\I - Zjv)), 

which is well-defined for Z = {Zi, . . . , Zjy) G Bi{8) since Z^ must be a strict con- 
traction < 1) for every k = 1, . . . ,N . In addition, C establishes a one-to-one 
correspondence from Bi{8) onto H^{8). The Szego kernel on Bi{8) can be transported 
on Hj^{8) by the Cayley transform. Thus, we introduce the Szego kernel on H^{8) by 
the formula: 

S{W, W) = F;,Fw', W, W G H48), 

where Fw = 2dia.g{{—i+W^))Sc-^(w)- Much more remains to be done in this direction. 
For instance, some classes of orthogonal polynomials of Jacobi type and their generating 
functions are considered in [7j. 

Finally, we mention that there are examples of polynomial relations for which there 
are no orthogonal polynomials. Thus, consider 

A = {X^-Xk\k = l,...,N}U {XkXi + XiXk I fc, / = 1, . . . , 2iV}, 

then TZ{A) ^ A(C^), the exterior algebra over C^. If is a unital q-positive definite 
functional on 71{A), then = for = 1, . . . , iV. This and the q-positivity of (p 

force (p{Xa-) = 0, therefore there is only one q-positive functional on 71{A) which is not 
strictly q-positive. Therefore there is no theory of orthogonal polynomials over this A. 
However, the situation is different for A = {X^Xi + XiX^ \ k,l = 1, . . . , 2A^}. This 
and other polynomial relations will be analysed elsewhere. 
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